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ABSTRACT 

Massive satellite accretions onto early galactic disks can lead to the deposition of dark matter in 
disk-like configurations that co-rotate with the galaxy. This phenomenon has potentially dramatic 
consequences for dark matter detection experiments. We utilize focused, high-resolution simulations 
of accretion events onto disks designed to be Galaxy analogues, and compare the resultant disks to 
the morphological and kinematic properties of the Milky Way's thick disk in order to bracket the 
range of co-rotating accreted dark matter. We find that the Milky Way's merger history must have 
been unusually quiescent compared to median ACDM expectations and therefore its dark disk must be 
relatively small: the fraction of accreted dark disk material near the Sun is about 20% of the host halo 
density or smaller and the co-rotating dark matter fraction near the Sun, defined as particles moving 
with a rotational velocity lag less than 50 km/s, is enhanced by about 30% or less compared to a 
standard halo model. Such a dark disk could contribute dominantly to the low energy (of order keV for 
a dark matter particle with mass 100 GeV) nuclear recoil event rate of direct dectection experiments, 
but it will not change the likelihood of detection significantly. These dark disks provide testable 
predictions of weakly-interacting massive particle dark matter models and should be considered in 
detailed comparisons to experimental data. Our findings suggest that the dark disk of the Milky Way 
may provide a detectable signal for indirect detection experiments, contributing up to about 25% of 
the dark matter self-annihilation signal in the direction of the center of the Galaxy, lending the signal 
a noticeably oblate morphology. 

Subject headings: Cosmology: theory — galaxies: formation — galaxies: evolution 



1. INTRODUCTION 

Prospects for the direct detection of dark mat- 
ter depend crucially on the phase space d istribution 
of dark matter near the Sun (Smith & L ewinI 119901 : 
lJungman et"all 119961: IVogelsberger et all |200^~^Unfor- 
tunately, detailed predictions for the dark matter distri- 
bution around the Earth are extremely difficult to con- 
struct from first principles; they require an understand- 
ing of spatial clumping on the scale o f the Solar S ystem 
(jKamionk owski fc KoushiappasI [20081 : [Peter '2009). and 
are almost certainly affected by poorly-understood bary- 
onic proces ses like the f ormation of the Galactic disk. 

Recently. [Read et all (|2008ll2009( l emphasized that the 
process of cosmological disk galaxy formation can sig- 
nificantly alter the dark matter distribution compared 
to canonical predictions that rely on simulations of the 
standard halo model. During the process of hierarchi- 
cal structure formation, merging satellite galaxies can 
get dragged into the plane of their host disk and deposit 
their dark matter in a structure, dubbed the dark disk 
([Read et al.l l2008f ). that is co-rotating with the Milky 
Way stellar disk and morphologically resembles a thick 
disk. If the ACDM cosmology represents the correct 
model of structure formation in the universe, it is certain 
that dark disks are virtually ubiquitous in disk galaxies. 
However, this contribution relative to the smoother halo 
component will depend sensitively on the formation pro- 
cess of each galaxy individually. In this work, we pro- 
vide the first focused attempt at constraining the dark 
disk contribution in the Galaxy, improving upon the ini- 
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tial di scussion of the Milky Way dark disk in [Read et al.l 
(2008) by comparing our simulations to the detailed kine- 
matic and morphological properties of the Galactic thick 
disk. 

If a substantial dark disk exists within the Milky Way, 
then it has dramatic implications for the direct and in- 
direct detection of dark matter, should that species con- 
sist of Weakly Interacting Massive Particles (WIMPs), 
which are pervasive in models of new physics at the weak 
scale. A co-rotating disk of dark matter will leave tempo- 
ral modulation signals in terrestrial nuclear-recoil exper- 
iments. It will also increase WIMP capture in the Sun 
and Earth and enhance the resultant flux of neutrinos 
from WIMP self-annihilatio n; both of th e se resu lts have 
been discussed in detail bv iBruch et al.l ()2009al |bh. We 
show in this work that a disk-like structure could be dis- 
cernible in the self-annihilation signal (depending on the 
total flux) from the center of the galaxy. Since the scale 
height of the dark matter disk will be larger than the thin 
disk scale height and smaller than any background halo 
flattening, such a morphological feature could provide an 
important handle on this indirect detection signal. 

The fraction of dark matter locked up in a co-rotating 
compon ent is expected to depend sensitively on merger 
history ([Read et al.ll2009( l. and this creates an important 
link between dark matter detection experiments and ef- 
forts in Galactic astronomy to constrain the accretion 
history of the Milky Way. Evidence is mounting that 
the merger history of the Galaxy is unusu ally quiescen t 
compared to typical ACDM expectations ([Wvsd [2009() . 
High-resolution simulations have shown that 1:10 mass- 
ratio accretion events, despite being fairly common in the 
ACDM paradigm (Stewart et al. 200 S), are irreco ncilable 
with the cold and thin Galactic disk (jPurcell et al.ll2009a[ . 
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Fig. 1. — Morphological and dynamical properties of the final 
stellar disks at the solar neighborhood for each fiducial orbital infall 
condition: the disk scale height z^i in kpc {upper panel), and the 
total velocity dispersion crtot in km/s [lower panel) . The shaded re- 
gions indicate observational boundaries: the Galactic scale heights 
and their e rrors are drawn fr om the SDSS tomographic results 
reported bv lJuric et al.l II2008I) . while the dynamical tem perature 
ranges represent the thin-disk spread observed by Nordstro m et ahl 
II2004I ) for total velocity dispersion as a function of stellar age in the 
solar neighborhood (the region represents stars with ages between 
1 Gyr and 8 Gyr and the dotted line marks the velocity dispersion 
for stars with a median age of ~ 2 — 3 Gyr), and the thick-disk 
range corresponds to the dispersions reported by .Soubiran et al.l 
1(2003). The infall model denoted by 0°s represents the accretion 
event involving a slower satellite galaxy; the three 60° infall models 
are marked with prograde, retrograde, and gas, according to the 
methods outlined in Sl2]and as listed in Table [Tl 

hereafter PK B09). This conclusion is consistent with 
the re sults of Read et al.l (|2008 ) and , Villalobo s & Helinl 
(|2008f ). who used simulations to show that 1:10 merg- 
ers generate heated systems that are grossly consistent 
galactic thick disks, but that are clearly incommensu- 
rable with the dominant thin disk o f the Galax y. More 
recently, the simulations of Mostcr et al.l (120090 showed 
that while these common mergers may be consistent with 
the broad population of galaxies, they produce remnants 
that are clearly thicker than the Milky Way within ~ 500 
pc of the disk plane. Empirically, the Galaxy appears 
to be deficient in stellar mass and angular momentum 
when balanced again st a sample of local spiral galaxies 
(jHammer et al.l[2007( ). further motivating the need for a 
focused program that constrains the Milky Way's dark 
disk specifically, rather than relying on general cosmo- 
logical expectations. 

We extend the methodology of PKB09 and use high- 
resolution simulations to study dark disk production. 
Our main observational constraint comes from observed 
kinematic and spatial properties of the old (~ 10 Gyr) 
thick disk of the Milky Way, which compri ses at most 
^ 20 % of the total mass of the Galaxy (e.g. \Juric et al\ 
\200^ ). We initialize stellar disks that are consistent with 
the observed properties of the Milky Way today and con- 
sider the impact of fairly massive, cosmologically com- 
mon accretion events. This approach is conservative: by 
initializing a disk that is as massive as the Milky Way 
disk today, we are exploring a case that is both more 
resistant to heating and more efficient in capturing dark 



matter from the accreted satellite than any realistic pro- 
genitor of the Milky Way's thick disk. More generally, 
by focusing on the oldest, thickest, and hottest stellar 
component of the Galactic disk, we are able to provide a 
most conservative constraint on the Milky Way's merger 
history. 

The heating of the initial disk and the creation of the 
dark disk are intimately connected, since heating requires 
the transfer of kinetic energy from the incoming satellite 
to the disk. For a fixed set of orbital parameters, more 
massive accretion events produce more stellar heating, 
and they also deposit more dark matter into a disk-like 
configuration. The transfer of kinetic energy depends on 
how close the satellite gets to the disk and its orbit about 
the disk; a closer and slower orbit results in larger heating 
and concomitantly larger stripping of dark matter from 
the satellite by the tides of the Milky Way. Thus, by 
comparing the heated stellar disk to the observed Milky 
Way thick stellar disk, we are able to bound the amount 
of coherently rotating dark matter that is deposited dur- 
ing the disruption of a large satellite galaxy. 

As alluded to above, the Milky Way has both a thin 
and thick disk (jGilmore fc Reidf ll983l ). with the latter 
component being significantly hotter and thicker than 
the former. By comparing the stellar remnants in our 
simulations with those of the Galactic thick disk implic- 
itly, we allow for the possibility that the thin disk of 
the Milky Way is regrown later through fresh accretion. 
While the regrowth of a thin disk could potentially act to 
make the heated disk thinner, it cannot make the heated 
disk colder. Thus our constraints based on the resultant 
velocity ellipsoid are the most robust. 

The majority of the accretion events we consider have 
mass ratios Msat:-/Whost = 1:10. This type of merger 
represents the dominant mode b y which dark mat- 
ter halos grow in ACDM models (|Purcell et al.l [20071 : 
IStewart et aT1l2008f l. They are also the most relevant for 
the formation of co-rotating dark disks, because more 
massive mergers are both quite rare and catastrophi- 
cally destructive to primary disks, while smaller merg- 
ers are less heavily affected by dynamical friction and 
leave behind much less dark matter (we explore such 
a case below). Thus the 1:10 events we consider here 
are both cosmologically common and capable of produc- 
ing a co-rotating dark disk without completely destroy- 
ing the primary galaxy. Note that while a number of 
less mass ive mergers may i ncrementally build up a thick 
disk (see lRead et al.l[2008| ). they will not produce a co- 
herently rotating dark disk unless the angular momenta 
contributed by the individual accretion events also co- 
incide constructively. This tendency toward a series of 
prograde accretions is not found in ACDM halos, at least 
for the few merger histories th at have been explored in 
detail (jKazantzidis et al.ll2008[ ). 

In addition to coUisionless experiments, we perform a 
hydrodynamical simulation in which the primary galaxy 
also hosts a gaseous disk, (a much more thorough 
treatment of these and similar results can be found in 
Kaz antzidis, Purcell, fc Bullock 2009, in preparation; see 
also lMoster et al]|2009f ) to determine whether this com- 
ponent plays a substantial role in the formation of a dark 
disk. We find that our results do not depend sensitively 
on the presence of gas as we have modeled it, though we 
caution that the treatment of ISM gas physics in galaxy 
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Fig. 2. — The three-dimensional dark matter density (top panels) and stellar density (middle panels) contributed by the accreted satellite 
galaxy for a selection of orbital conditions, viewed edge-on with respect to the primary galactic stellar disk and calculated in a vertical 
slice with a thickness of 1 kpc which passes through the center of the galaxy. For reference, the host galaxy's structure is mapped in the 
left panel of each figure group and the density maps of the heated primary disk stars are shown in the lower panels; each box is 20 kpc on 
a side and the solar position is marked with 0. 



simulations remains uncertain. Based on our implemen- 
tation, the presence of gas does not curtail the heating 
of a primary disk that occurs during dark disk creation. 

2. METHODS 

In a suite of high-resolution coUisionless simulations, 
PKB09 investigated the response of a galactic disk to the 
infall of cosmologically-common accretion events involv- 
ing satellites with one-tenth the mass of the host halo. 
For our primary system, we focus here on the Galaxy 1 
model from PKB09, which is a good match to the Milky 
Way today, and should provide a higher probability of 
dark disk creation via satellite dragging than any pre- 
sumably less-massive progenitor of the Galaxy at higher 
redshift. We now provide a brief overview of our simula- 
tions and refer the reader to PKB09 for a more complete 
discussion of the methods and results regarding the mor- 
phological and dynamical changes undergone by stellar 
disks in response to these accretion events. 

Our primary galaxy is constructed according to the 
fully self-consisten t distribution functions prescribed by 
iWidrow et al.l ('2008) and is therefore an equilibrium so- 
lution to the coupled coUisionless Boltzmann and Pois- 
son equations. The model's stellar disk is initialized 



with an exponential scale length Rd = 2.84 kpc and a 
vertical distribution described by a sech^ function with 
scale height = 0.43 kpc, containing 10^ particles with 
a total disk mass A/disk = 3.6 x IO^^A/q. The mas- 
sive central bulge follows a Sersic profile of effective ra- 
dius Re = 0.58 kpc and index n = 1.118, and has a 
mass Afbuigc = 9.5 x IO^A/q distributed among 5 x 10^ 
particles. These stellar components are embedded in a 
dark host halo composed of 4 x 10® parti cles which fol- 
low th e canonical NFW density profile of iNavarro et al.l 
(|1996f ). with scale radius Vs = 14.4 kpc and virial mass 
AThost — lO^^Af0. This particular set of parameters was 
chosen in order to minimize secular effects such as bar 
formation, as well as artificial heating induced by the 
interaction of disk particles with more massive halo par- 
ticles. 

The abscissa of Figure [T] as well as the top row of 
Table [T] denote the infall models we consider. Every 
simulation tracks the infall of a satellite galaxy mod- 
eled within a subhalo of virial mass Afgat — 10^^ Mq at 
z = 0.5. The satellites are each composed of 9 x 10^ 
particles in a density structure well fitted by an NFW 
profile with concentration Cvir — 14, and contain a stel- 
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Fig. 3. — Rotational velocity histograms at the solar radius for accreted dark matter (thick solid), accreted stars (shaded region), and 
host halo dark matter (thin dotted) for infall inclinations corresponding to those shown in Figure [31 calculated with respect to the local 
standard of rest defined by the rotation speed vi^sR of the primary galactic stellar disk. The dark matter values are normalized in each bin 
to the total number of dark particles in the area defining the solar neighborhood; the number of accreted stars in each bin is normalized to 
the total number of stars in the same region. In each panel, the thin solid line represents the sum of the two Gaussians representing dark 
matter in the host halo and from the accreted subhalo. 



lar mass — 2.2 x IO^Mq distributed among 10'"' par- 
ticles in a spheroid with Sersic index n ^ 0.5. The stel- 
lar mass es for the satell ite galaxies are drawn from the 
model of lConrov fc Wechsler, (,2009l) . which describes the 
redshift-dependent connection between the observed spa- 
tial abundance of galaxies and the corresponding abun- 
dance of predicted dark matter halos. Fiducially, we 
investigate five orbital infall inclinations, including four 
prograde orbits with 9 = (0°, 30°, 60°p, and 90°) and one 
retrograde 60°r orbit. The satellites are initialized 120 
kpc from the host halo center, and we set our initial or- 
bital vectors according to the distribution of subst ructure 
accretions drawn fro m cosmological simulations ( BensonI 
l2005t[Khochfar &: Burkert 2006). with radial and tangen- 
tial velocity components equal to Vr = 116 km s~^ and 
vt = 77 km s~^ respectively. Additionally, in a model 
we designate 0°s, we investigate a case where the same 
satellite galaxy is traveling more slowly than in the fidu- 
cial 0° case, with Vr — 58 km s^^ and vt — 38.5 km ss^^; 
these orbit al parameters c orrespond roughly to the lower 
Icr limit of lBensonI (|2005D . 

To test the effect of a full treatment of hydrodynamics 
on the morphological and dynamical changes induced by 
an accretion event, we repeat the simulation involving a 
massive satellite galaxy infalling along a prograde orbit 
with an inclination of 60° (hence the designation 60°g), 
having converted a modest fraction fg — 15% of the ini- 
tial primary galaxy's star particles into a gaseous compo- 
nent with the same density structure as the stellar disk. 
Our hydrodynamical prescription includes atomic cool- 
ing for a primordial mixture of hydrogen and heliu m, and 
the sta r formation algorithm is based on the work of lKatd 
()1992|) . in which gas particles in cold and dense regions 
form star particles at a rate proportional to the local dy- 
namical time; star formation occurs when the gas density 
exceeds 0.1 cm~^ and the gas temperature drops below 
1.5 X 10"* K. Supernova feedback is imple mented accord- 
ing to the blast-wave model described in iStinson et al.l 
()2006[ ). in which the energy deposited by a Type-II su- 
pernova into the surrounding gas is 4 x 10^'^ ergs. This 
parameter set p roduces realistic galax ies in cosmologi- 
cal simulations (|Governato et al.ll2007t ). Again we refer 
the reader to Kazantzidis, Purcell & Bullock (2009, in 



preparation) for more details regarding the suite of ex- 
periments to which this current test case belongs. 

As a test of the fractional contribution made by a much 
smaller satellite galaxy, we simulate an additional ac- 
cretion event with identical orbital parameters to the 
fiducial infall, involving a subhalo of mass Mgat 4 x 
IO^^A/q, i.e. with mass ratio 1:25. For this case, we 
choose a planar subhalo orbit [0 — in order to maxi- 
mize the likelihood of dark disk formation. However, this 
experiment resulted in a trivial increase (< 1%) in the 
fraction of dark matter in the solar neighborhood, and is 
thus not listed in Table [T] 

3. RESULTS 

3.1. Gross Comparisons to the Galactic Stellar Disk 

We summarize in Figure [T] the observable properties 
of our initial and remnant stellar disks, with each model 
indicated by name along the horizontal axis. The upper 
panel shows the disk scale height z^, derived from fitting 
the minor axis surface density profile at a defined solar 
radius {Rq = 8 kpc), using the two-component form 

Y.{z) = i;dsech^(z/zd) + Sdiffusc sech^(z/zdiffuso)- (1) 

This decomposition provides a minimum estimate of the 
resultant disk scale height (z^) by allowing for a sec- 
ondary, much thicker component (zdiffusc), which domi- 
nates at large height and low surface brightness. Gen- 
erally, we find significant variance in these values with 
orbital inclination angle; the disk scale height Zd ~ 
1 — 2 kpc, while the faint component typically has 
•Zdiffusc ~ 4 — 7 kpc. Note that the initial thin disk of the 
primary galaxy has been destroyed: our minimal scale 
heights are comparable to or larger than the Galactic 
thick disk scale height, as demonstrated in the compar- 
ison of Figure (TJ where observed scale heights for the 
Milky Way thin and thick disks are shown as shaded hor- 
izontal bands. These values are drawn from f juric et al.1 
(2008), in which an exponential disk scale height is pre- 
sented; we convert this value here to a sech^ scale height 
such that the two vertical surface density profiles de- 
crease by roughly the same amount within 1 scale 
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Fig. 4. — Correlations between the kinematic properties of dark 
matter and stars accreted onto the primary galaxy during satellite 
infall. The upper left panel shows the amount Vi^^ by which each 
species of accreted material lags the rotational speed of the primary 
galactic disk; negative (positive) quantities indicate faster (slower) 
rotation than the local standard of rest as defined by the mean stel- 
lar velocity of disk stars in the solar neighborhood. The remaining 
three panels plot each axis of the velocity ellipsoid for both species. 
In all figures, the ordinate marks the values for accreted stars, 
while the abscissa denotes the corresponding values for accreted 
dark matter; the dotted line defines where the two values are equal, 
while the solid line forms the upper bound a^^^ = 1.25(7''''""'' of 
an envelope which roughly contains the simulation results. 

height of the disk planc^. The total three-dimensional 
velocity dispersion atot — + <^e + ^"^Y^"^ of each disk 
at the solar neighborhood is presented in the bottom 
panel of Figure [TJ The significant increase in this quan- 
tity from the initial model to the resultant disk betrays 
the extraordinary heating undergone by the primary disk 
during the merger. The shaded bands again reflect ob- 
served velocity dispersions for the Milky Way thin and 
thick disks, as drawn from iNordstrom et al.l ()2004[ ) and 
iSoubiran et al.l ()2003D . respectively. 

If we associate our remnant disks with the Milky Way 
thick disk, most of our resultant systems are thicker and 
hotter than the Galaxy. This comparison is particularly 
constraining for two reasons. First, the age distribution 
of stars in the thin disk includes stars that are 8-10 Gyr 
old, suggesting than any merger of this kind must have 
happened well before z ~ 1 (when the primary disk was 
even smaller and less able to cause significant dragging 
into the disk plane). Moreover, while regrowth of a new 
disk could potentially reduce the scale height of the thick- 
ened disk, it would not reduce the velocity dispersion 
of the stars. Only the prograde and retrograde 60° or- 
bits and the 90° degree orbits result in a system that 

^ Note th a t mu ch of the relevant literature, including e.g. 
IMoster et all ||2009| '). uses the typical conversion factor of two, 
which is only valid at heights much larger than the scale height 
of the disk. Our conversion method captures more accurately the 
behavior of the surface brightness profile at heights comparable to 
and smaller than the disk scale height. 
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Fig. 5. — Similar to Figure [T] morphological and dynamical 
properties of the final accreted dark disks at the solar neighbor- 
hood for each fiducial orbital infall condition: the local density 
ratio Pdl Ph (upper panel), and the fraction of local dark matter 
moving slowly (|fiag| < 50 km/s) with respect to the stellar disk's 
rotational speed (lower panel). In the latter case, the thin dotted 
histogram denotes the contribution to the fraction of slow-moving 
dark matter made by the initial host halo, while the colored points 
indicate the sum of this component and the amount of co-rotating 
dark matter deposited by the accreted satellite galaxy. The shaded 
regions correspond to the orbital infall models that produce stellar 
remnants with morphology and kinematics similar to the Galactic 
thick disk, as demonstrated by Figure [T] 

is marginally viable, both morphologically and dynami- 
cally. 

3.2. Morphological and Kinematic Characterizations of 
the Dark Disk 

Each of the merger simulations provides an accreted 
dark component and an accreted stellar component, but 
the nature of these components depends sensitively on 
the interaction, as illustrated in Figure [D for three of our 
simulations. For reference, the upper left panel presents 
a central slice of the dark matter density in the pri- 
mary galaxy's host halo prior to the accretion, the mid- 
dle left panel shows a slice in the stellar distribution in 
the primary galaxy, and the lower panels represent the 
heated stars belonging to the original disk only. The 
three columns on the right for the upper two rows display 
the resultant distributions of accreted dark matter and 
accreted stars; as might be expected, the lower-latitude 
event produces the most disk-like accreted dark matter 
morphology, with a qualitatively similar accreted stel- 
lar distribution that is somewhat more rotationally sup- 
ported owing to the fact that the stars are more tightly 
bound in the satellites than the dark matter. As noted 
in Table [U the low-latitude simulations also produce the 
largest fractional density pa/ Ph of accreted dark mat- 
ter compared to background halo dark matter in the 
solar neighborhood^. High-latitude events create less 
discernible dark disks at the solar position and deposit 
proportionally less dark matter there; we also note that 

Throughout this work, we define the solar neighborhood as 
the area between galactocentric radii 7 < i? < 8 kpc and bounded 
vertically by |z| < 2 kpc; we also note that variance in these bound- 
aries does not lead to significantly altered results. 
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our test model involving the planar infall of a satellite 
galaxy with mass ratio Mgat^Afhost = 1:25 evolves over a 
much longer timescale than that of the fiducial models, 
and also results in a negligible contribution to the local 
dark matter quotient, indicating that even a multitude 
of small subhalos with similar orbits may not be able to 
form a significant dark disk. For the models with mass 
ratio Msat:Afhost = 1:10, Table [T] provides dark disk and 
accreted stellar disk scale heights at the solar neighbor- 
hood, as computed by fitting a sech^ profile. 

An additional characterization of interest is the veloc- 
ity distribution of the dark disk component. It is useful to 
parametrize the dark matter in the solar neighborhood 
of our simulations with a double Gaussian distribution 
that represents the original host and accreted subhalo 
material: 

fivR.V^.V:,) = /host(wi?,W0,Wz) +/acc(Wfl,1'0,Wz), (2) 

where in each term 

f{vR,v^.v,)=Ae-^^^^-'^^^"'^^l (3) 

The sum involves each velocity component, i = R^cji^z. 
We work in a frame that co-rotates with the resultant 
stellar disk such that the host halo produces a net dark 
matter wind in the direction that is close to the pri- 
mary galaxy rotation speed (w0)host ^ —200 km/s and 
v'r = Vz '-^ 0. Accreted matter can co-rotate or anti- 
rotate with the stellar disk and we quantify this rotation 
in accreted material by the lag speed wiag = —v^, such 
that a zero lag corresponds to precise co-rotation, while 
a positive lag trails the stellar disk and a negative lag 
means that the dark disk is rotating faster than the stel- 
lar disk. The resultant velocity distribution parameters 
for each accreted component in our simulations are pre- 
sented in Table [TJ 

In Figure [31 we show the collapsed distributions in ro- 
tational velocity for the background host halo as well 
as accreted dark matter, along with our double-Gaussian 
fit to the total amount of dark matter in the solar neigh- 
borhood, for the same three resultant galaxies that were 
illustrated in Figured) 30°, 60°p, and 60°g. The cen- 
ter panel also presents the velocity distribution for the 
60°r retrograde merger. The filled histograms in each 
panel show the velocity distributions for accreted stars, 
normalized relative to the total stellar mass at the so- 
lar location. Two points are immediately apparent; first, 
that lower-latitude accretion events produce smaller lag 
velocities in both accreted species, and secondly that the 
lag speed of accreted stars grossly mimics that of the ac- 
creted dark matter for accretion events at all latitudes. 

3.3. Correlating the Dark and Stellar Disks 

The kinematic relationship between the accreted dark 
and stellar disks (a phenomenon also reported by 
iRead et al]|2008l |2009( ) is important because it may pro- 
vide an avenue to constrain the dark disk velocity dis- 
tribution via observational studies that can isolate an 
accreted stellar component. In Figure SJ we show the 
correlations between accreted dark matter and accreted 
stars in each velocity dispersion component, as well as 
the rotational lag speed, which correlates very well be- 
tween the two accreted species. We also note that the 
velocity dispersion of accreted stars generally provides a 
lower limit on the velocity dispersion of the dark disk in 



each directional component; the correlations in Figure [4] 
indicate that dark matter dispersions are consistently at 
least a factor of 1.25 larger than their stellar counter- 
parts. 

Though not explicitly shown here, we also find a con- 
nection between the velocity-ellipsoid axis ratio Uz/cfr of 
the accreted dark matter and that of the stars formerly 
belonging to the satellite galaxy, in that the two collision- 
less events involving the high inclination angle of 6* = 60° 
both have vertical to radial velocity dispersion ratios of 
approximately unity for all accreted material, while the 
velocity ellipsoids resulting from the low-latitude events 
are significantly more oblate, with both stellar and dark 
axis ratios roughly equal to ~ 0.5: a similar value to that 
obtained obscrvationally for the Galactic so lar neigh- 
borhood (,Gomez et _aL. 1990i: iNordstrom et al". 2004). It 
should be noted here that the polar infall results in a 
final velocity-ellipsoid axis ratio coincident with that of 
the low-latitude events, due to a very large dispersion in 
the radial velocity of both accreted species in this case. 

3.4. Summary of Results 

We use two simple metrics to characterize our main re- 
sults in Figure O as in Figure [TJ each model is indicated 
along the horizontal axis. The upper panel presents the 
accreted dark matter fraction Pd/ Ph in the solar neigh- 
borhood; the lower panel presents the fraction of local 
dark matter that is roughly co-rotating with the stellar 
disk, which is potentially the most salient characteris- 
tic of an accreted dark disk. Specifically, we define the 
latter quantity to be the fraction of accreted dark parti- 
cles with rotational speeds within 50 km/s of the stellar 
disk's velocity, i.e. \Vfi, — vlsr\ < 50 km/s. Note that 
by this definition, even the initial spherical halo has a 
non-zero co-rotating dark matter fraction; for reference, 
the dotted line in the lower panel of Figure [5] shows the 
fraction of dark matter particles initially from the back- 
ground halo that are co-rotating with the stellar disk in 
the final simulation snapshot. 

As might be expected, the most significant dark disks 
form during the planar accretion events 0° and 0°s, which 
respectively produce 32% and 35% accreted dark mat- 
ter fractions in the solar vicinity and factors of 1.7 and 
1.4 enhancements in the co- rotating dark matter frac- 
tion. However, these planar events also produce stellar 
disks that are much hotter than the thick disk of the 
Milky Way (as shown in Figure [1]), indicating that these 
cases are not appropriate baseline models for the Galaxy. 
Though not evident in Figure (TJ the slower orbit actually 
produces more radial heating than the fast orbit (see Ta- 
bic [T]) and is therefore more discrepant with the Milky 
Way thick disk velocity structure than the standard or- 
bit. 

The three shaded bands in Figure [5] highlight the three 
model cases that produce stellar disks marginally consis- 
tent with the thick disk of the Milky Way, when con- 
sidering morphology as well as dynamical temperature. 
The most optimistically viable model for a dark disk is 
the 60-degree prograde orbit, which results in a 16% ac- 
creted dark matter contribution and a 14% co-rotating 
dark matter fraction; the latter quantity represents a 
~ 30% increase compared to the co-rotating dark matter 
fraction in the smooth halo case. As noted in Table [1] 
the dark disk scale height in this case is Zd = 4.6 kpc and 
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lags the LSR with viag = 46 km/s. 

Although the thick disk we produce is obviously more 
massive than the Galactic thick disk, the similarity in 
morphology and kinematics motivate us to move beyond 
limits, and speculate on the properties of the actual dark 
disk in the Milky Way. Given the r ange of parameters 
probed here and in lRead et al.l ()2008[ ). and operating un- 
der the assumption that the Milky Way's thick disk was 
formed in an event similar to the type we simulate here, 
we expect that the real Galactic dark disk contributes 
roughly 10% to the local dark matter density. As we dis- 
cuss below in fJSj such a dark disk does not significantly 
enhance the possibility of direct detection of dark matter; 
however, it will be important to include in any detailed 
analysis aimed at interpreting an observed WIMP detec- 
tion signal. Moreover, such a dark disk may provide a 
detectable morphological signature towards the Galactic 
center in experiments designed to indirectly detect dark 
matter via its self-annihilation products. 

4. COMPARISON TO PREVIOUS WORK 

In a re cent effort t o quan tify disk heating in a ACDM 
context, iRead et ahl (|2008f) used simulations similar to 
those presented here in order to show that massive satel- 
lite accretion events involving a Galactic-type thin disk 
not only produce thick stellar systems but also result 
in the deposition of accreted dark matter into a disk- 
like component that co-rotates to some degree with the 
primary galaxy. Specifically, they analyzed dark disk 
formation during the infall of various subhalo types (of 
particular interest here are their LMC and LLMC mod- 
els, representing satellite galaxies with total mass Mgat — 
2.4 X 10^° and 1.0 x 10^^ respectively), finding that for 
low-inclination mergers, the resulting density ratio of ac- 
creted dark matter to that of the background host halo in 
the solar neighborhood was Pd/ Ph = 0.22 for the smaller 
subhalo and 0.42 for the more massive satellite that cor- 
responds closely to our fiducial infalling system. In the 
follow-up analysis of cosmological simulations of Milky- 
Way- type disk galaxies (|Read et al.li2009i) . the authors 
found similar results for systems with a wide range in 
accretion history activity, showing that several massive 
mergers at late times can result in a dark disk with local 
density similar to that of the host halo itself. 

Although our simulated phenomena are qualitatively 
similar to those presented by Read ct al. (2008), there 
are two key points of distinction to be drawn. The two 
methods are only significantly discrepant in the choice 
of initial radii and velocities of the infalling satellites. 
Our work adopts satellite galaxy velocity vectors drawn 
from subhalo infall distributions found in cosmological 
simulations, and thus we set the initial position of the 
subhalo far (~ 120 kpc) from the host halo's center in 
order to minimize the sudden change in potential felt by 
the primary galactic disk and also to imitate the cos- 
mological velocity conditions, which are measured at the 
host's virial radius. In contrast, the infalling satellites of 
iRead et aD ([20 OS) begin much closer (~ 30 kpc) to the 
primary disk's center and involve subhalos traveling at 
speeds slower than our fiducial initial velocity by roughly 
fifty percent; slightly faster than our 0°s case. The au- 
thors motivate the choice of a slow orbit at a small initial 
Galactocentric apocenter by suggesting it replicates the 
orbital parameters of satellites found in a loose group 



environment that is accreted onto a Galaxy-sized host. 
However, it still remains to be shown that such a velocity 
vector at ~ 30 kpc can arise naturally for a massive satel- 
lite within a cosmological setting. Moreover, any process 
that would act to bring a subhalo to such a low energy 
state in the disk plane would almost certainly heat the 
primary disk to a degree unallowable by observations. 

The difference in initial conditions may be largely re- 
sponsible for the small systematic differen ces between 
our simulated dark disk results and those of iRead et al.l 
([2(308). While our maximal local density ratio is Pd/ Ph ~ 
0.35, they find 0.42 in their experiment labeled LLMC- 
10°, an accretion event involving a subhalo of similar 
mass. The methodological variance may also explain why 
our disks are slightly dynamically hotter. Our 0° case has 
ctot ~ 117 km/s and our 0°s case has crtot ^ HO km/s, 
while their LLMC-10° model has dtot ~ 105 km/s. It is 
important to emphasize, however, that both sets of sim- 
ulations find resultant thick disks that are hotter than 
the observed thick disk of the Milky Way, which has a 
total dispersion of just ^ 85 km/s. Therefore, in a broad 
sense, we agree that these disks are not ideal analogues 
to the Milky Way. 

We note that a more substantial dark disk component 
could emerge following multiple similar accretion events, 
as investigated in cosmological s imulations of Galaxy- 
sized disks by IRead et all (|2009f ). but such a series of 
late-time mergers cannot be reconciled with the rela- 
tively thin disk of the Milky Way. Unfortunately for 
experiments aimed at detecting local dark matter, the 
quiescent Galactic accretion history favored by our re- 
sult indicate that the relevant observable quantities are 
probably far less affected by a dark disk than we might 
have hoped. 

5. DISCUSSION AND INTERPRETATION 

In our experiments, we have adopted assumptions 
throughout that maximize the likelihood for producing a 
dark disk remnant and for preserving a thin, cold disk by 
initializing a primary disk that is as massive as the Milky 
Way disk today. This approach is quite conservative be- 
cause we compare the heated stellar remnant to the old 
thick disk of the Milky Way, which is at leas t five times 
less m assive than our primary system {e.g. \Juric et al\ 
\200A . and references therein). We have argued that the 
dark disk formed in our 60° prograde merger case pro- 
vides a fair limit on the properties of an underlying ac- 
creted dark component of our Galaxy, with a local con- 
tribution Pd/ Ph — 0.15 that is at the low end o f dark disk 
contributions estimated bv ' Read et al.l (|2008f ) for Milky 
Way type galaxies in ACDM. 

Based on this limit, we expect that nuclear recoil ex- 
periments designed to directly detect local dark mat- 
ter would only receive small boosts to the event rate, 
preferentially at lower energies. Dark disk fra c tions o f 
Pd/ Ph ^ 10% were explored by iBruch et all (|2009af) . 
who found an order unity enhancement in the differential 
event rate at keV recoil energies for a 100 GeV dark mat- 
ter particle. For a TeV mass dark matter particle, such 
an enhancement would be present at higher energies (ac- 
cessible to CDMS-II) . While an order unity enhancement 
will no t significantly impact the event rate, iBruch et al.l 
((20093) showed that the phase of the annual modulation 
signal is sensitive to both the dark disk and the mass of 
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Fig. 6. — The dark matter annihilation signal in Galactic coordinates, mapped with a color logarithmically proportional to the integrated 
square of the particle density along a line of sight toward the center of the galaxy, following the prograde satellite galaxy infall with 
inclination 9 = 60°. Note that the combined signal {right panels) is noticeably more oblate than the host material alone {left panels; shown 
after the accretion event); the accreted dark matter {center panels) contributes roughly ~ 10 — 25% of the annihilation rate at Galactic 
longitudes ~ 5° — 10°. In the lower set of figures, we test the effect on the total annihilation signal for a hypothetical dark disk with a local 
density ratio at the solar neighborhood Pd/ph ~ 1; this boost is also an 0(1) effect, contributing ~ 50% of the total signal over the same 
range in longitude. We have argued that such a large dark disk component would likely have to form through secular processes unrelated 
to satellite accretion events. 

the dark matter particle. The phase of the annual mod- 
ulation signal produced by the dark disk depends on the 
motion of the Sun relative to the dark disk. The phase 
of the modulation of the total event rate in a given en- 
ergy window depends on both the mass of dark matter 
particle and the dark disk fraction. If the dark disk frac- 
tion is bracketed from the lower end, this will result in 
a prediction for the phase given a dark matter particle 
mass. This may be testable at future detectors if a pos- 
itive signal is seen. In addition, these results also imply 
that directional detectors should see the WIMP wind di- 
rection change as a function of recoil energy, the details 
of which would depend on the energy window, the mass 
of the WIMP, the dark disk fraction and the motion of 
the Sun with respect to the dark disk. 

An increase in the local number density of particles 
that co-rotate with the Sun will also enhance WIMP cap- 
ture in the Sun and the Earth. These captured particles 
self-annihilate into standard model particles, including 
neutrino pairs. Thus a larger co-rotating fraction of dark 
matter will result in a larger neutrino flux potentially 
observable with experiments like Super-Kamiokande and 
IceCube. In detail, the capture rate depends sensitively 
on pd/ Ph, the dark disk lag speed wiag, and the veloc- 
ity dispersion of the dark disk. The timescale for cap- 
ture of WIMPs at Earth is small compared to the age 
of the planet and hence the WIMP density in the Earth 
has not yet reached equilibrium. In the Sun, as shown 



by iBruch et al.l ()2009bf ). the situation is the opposite 
for most regions of parameter space and the number of 
WIMPs would have reached an equilibrium value. Thus 
the annihilation rate of WIMPs inside the Earth is pro- 
portional to the capture rate squared, while in the Sun 
it is equal to half the cap ture rate (given the equilib- 
rium condition there). Br uch et al.l ()2009bl ) show that 
this leads to two to three orders of magnitude enhance- 
ment in the flux from the Earth and an order of magni- 
tude enhancement in the flux from the Sun for pd/ Ph = 1 
and an assumed isotropic Gaussian velocity distribution 
for the dark disk with lag wiag = 50 km/s and 1-D disper- 
sion ad — flag- They also consider dark disk parameters 
Pd/ph — 0.25 and <7d — 100 km/s and state that this 
does not lead to a large boost in the signal, owing pri- 
marily to the large velocity dispersion. By comparison, 
our simulated prograde accretion event with 9 = 60° re- 
sults in a local density ratio of 0.15, velocity dispersion 
~ 90 km/s, and similar lag speed vug ~ 50 km/s. 
Thus we do not expect large boosts to the neutrino flux 
from the Earth or the Sun, and our results imply that 
the likelihood of WIMP detection through this indirect 
channel is not significantly increased due to an accreted 
dark disk. 

Although the prospects for direct detection and indi- 
rect detection via WIMP capture are not significantly al- 
tered by the accreted dark disk, we note that the accreted 
dark matter distribution is relatively more prominent to- 
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ward the center of the galaxy, where most of the satellite 
galaxy's material settles. The center of the Galaxy is 
a prime target for experiments aimed at detecting dark 
matter indirectly via h igh energy annihilation products 
(jBergstrom et al.l Il998l ). as the flux scales proportion- 
ally to the square of the WIMP density. In Figure [H 
we view the halo's center from the vantage point of the 
solar neighborhood, after the prograde 60° satellite in- 
fall, and estimate the annihilation signal at Earth by 
integrating along lines of sight. We notice immedi- 
ately that while the host's dark matter has retained ap- 
proximate sphericity (left), the combined signal (right) 
which includes the accreted subhalo material (middle) 
has a significantly more oblate contour. Approximately 

10 — 25% of the total contribution at galactic lon- 
gitudes between 5° and 10° comes from accreted dark 
matter. As in the case of direct detection, the dark disk 
does not significantly change the likelihood of detection; 
however, if there is a positive detection, then our results 
provide motivation to search for a disk-like component 
with scale height and length that are different from that 
of the Milky Way thin disk. 

Our focu s in th is work has been on the scenario of 
iRead et al] ()2008D . who showed that dark disks could 
be created via satellite accretion events. We have ar- 
gued that under this scenario, the relative importance of 
the accreted dark disk can be constrained via detailed 
comparison to the properties of thick disk stars in the 
Galaxy and that this constraint limits the dark disk's 
density contribution to ~ 0.2phaio in the solar neighbor- 
hood. However, it is interesting to consider the theoret- 
ical possibility that a dark disk with density ratio 0(1) 
could form in response to some other process that did 
not heat the disk. For example, one could imagine that 
large infalling gas clouds could transfer angular momen- 
tum to the dark matter, creating a disk-like dark compo- 
nent without disturbing the primary disk significantly. 
Such a process may not be easy to arrange, given the 
large amount of angular momentum transfer required to 
produce a co-rotating dark disk as well as the fact that 
observed disk galaxies seem to require something close 
to angular momentum conservation. 

Nonetheless, it is an instructive exercise to determine 
the effect of an 0(1) dark disk on indirect detection sig- 
nals. To investigate this, we have scaled the local density 
ratio for the 60° prograde case such that Pd/ Ph — 1, and 
the resultant annihilation signal's sky map is shown in 
the lower panels of Figure El The disk-like nature of 
the signal is clearly apparent and could be easily distin- 
guished if any dark matter self-annihilation is observed 
toward the center of the Galaxy. Therefore, if such a 
dark disk could somehow form without involving satel- 
lite accretion events, it would be discernible through in- 
direct detection experiments. Coupled with the boost in 
the direct detec tion signal that is e xpected for such an 
0(1) dark disk (|Bruch et al.ll2009b[ ). the jointly observ- 
able consequences of such a component would be signifi- 
cant. 

In a similar vein, it is worth noting that the quantity 
of interest for direct detection experiments is the fraction 
of slow-moving dark matter particles. From Figure[5l we 
see that the 30° case shows an order unity enhancement 
in slow-moving dark matter particles, and Figure [1] shows 
that this run produces a disk that is only marginally hot- 



ter than the thick disk of the Milky Way. However, the 
scale height of the disk is considerably larger. One may, 
however, hypothesize that a later stage of disk formation 
could change the disk scale height substantially without 
heating the stars further. Therefore, we urge the reader 
to keep in mind that 0(1) enhancement in slow-moving 
particles is possib le in the accreted dark disk scenario of 
iRead et al] (|2008D . but that the expected velocity dis- 
persion in the 4> direction is still large - in our 30° case, 
it is 83 km/s (see Tabled]). 

In this work we have argued th at within th e context of 
the accreted dark disk scenario of lRead et al.l HOOH), it is 
likely that the dark disk of the Milky Way contributes ap- 
proximately 10 — 20% to the local dark matter density. If 
so, then its presence may be important to include in any 
attempt to interpret a detection signal of WIMP dark 
matter in the solar neighborhood and for the indirect 
detection signal from the Galactic center. Given this, it 
will be important to con strain the d a rk dis k to a higher 
degree of accuracy. Both IRead et al.l (|2008f l and we have 
demonstrated that there is a significant dynamical rela- 
tionship between dark matter accreted during the infall 
of a massive subhalo and the stellar mass contributed 
by that satellite galaxy to the primary system (Table [1] 
and Figure |4]). Therefore one can anticipate constraining 
the properties of any Milky Way dark disk by isolating a 
dominant subpopulation of accreted Galactic disk stars, 
perhaps by some combination of chemical and dynami- 
cal tags. For example, if we were to consider the unlikely 
scenario where the entire thick disk of the Milky Way was 
deposited by a single accretion event (i.e. it is devoid of 
any pre-existing disk stars), then we may consider the dy- 
namical properties of the Galaxy's thick disk as a means 
to constrain the dark disk: (ctr, ff^, (Tz) = (63 ib 6, 39 ± 
4, 39 ±4) km/s, ?;iag = 51.5 km/s rSoubir an et al]l200l . 
and Zd = 0.90 ± 0.18 kpc (fjuric et al. 200i). Under this 
interpretation, our simulation results would imply that 
the lag speed of the MW dark disk is ^ 50 km/s and that 
it must be quite hot, with 1.25 (cta', cr^, cr^) ~ (88, 51, 62) 
km/s, and with a very thick > 2 kpc. These numbers 
are almost certainly lower limits; as we have shown in 
Table [U accreted stellar distributions tend to be both 
hotter and thicker than the primary disks they heat. If 
there is any pre-existing disk population associated with 
the thick disk of the Milky Way, the accreted portion 
is almost certainly hotter than the composite thick disk 
values used in our analysis. 

Although theoretical predictions based on galaxy for- 
mation models as they relate to dark matter detection 
are still in the nascent stage, it is clear that the stan- 
dard halo model is insufficiently equipped to allow precise 
predictions for measurable quantities of interest; scat- 
tering rates and flux scalings depend crucially on the 
presence of locally coherent substructural flows. Obser- 
vations have thus far been unable to accurately estimate 
the degree to which the Galactic thick disk's population 
is composed of stars stripped from a satellite galaxy dur- 
ing tidal disruption, as opposed to material heated to 
ejection from the t hin stellar disk (as investigated by 
iPurcell et al.|[20b9bf ). In principle, however, a combina- 
tion of high-precision kinematic measurements and de- 
tailed chemical composition data should be able to de- 
fine a phase-space in which accreted and original stars 
occupy disparate regions. Placing such dynamical con- 
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straints on the population of stars obtained by the Milky 
Way during a significant accretion event may enable us 
to identify the kinematics of the dark matter deposited 
during that same event, thereby placing more stringent 
limits on the availability of a dark Galactic disk for the 
purpose of detection experiments. 
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TABLE 1 

Properties of accreted dark matter and stars in the solar neighborhood 



Quantity Orbital Inclination 



(DM) 


e = 0° 


Q = 0°-slow 


e = 30° 


e = 60° 


= 60° -gas 


e = 60°-retro 


e = 90° 


Units 


Pd/ph 


0.32 


0.35 


0.17 


0.16 


0.29 


0.21 


0.11 






132.7 


226.5 


150.5 


191.6 


157.4 


193.5 


483.0 


[km/s] 




78.2 


106.2 


83.3 


88.9 


87.1 


99.7 


75.1 


[km/s] 




63.4 


73.1 


78.1 


181.2 


152.3 


185.9 


241.2 


[km / s] 


"lag" 


-11.7 


24.1 


-9.6 


46.1 


51.5 


339.6 


202.2 


[km/s] 


Zd 


2.6 


1.8 


4.2 


10.9 


9.2 


9.2 


36.6 


[kpc] 


(stars)'' 




91.5 


119.9 


77.6 


127.6 


107.2 


143.7 


230.1 


[km/s] 




61.8 


57.0 


63.9 


53.6 


60.6 


61.4 


34.9 


[km/s] 




39.3 


35.7 


41.9 


138.9 


104.3 


196.1 


88.8 


[km/s] 


flag 


-9.8 


21.3 


-1.6 


68.3 


49.5 


286.2 


210.6 


[km / s] 


2d 


1.7 


0.8 


1.6 


4.6 


2.8 


9.4 


11.5 


[kpc] 



"Lag velocities are calculated in the frame of the local standard of rest (LSR), defined by the mean rotational speed of the primary galaxy's 
disk stars at the solar location; positive (negative) values indicate a slower (faster) rotation than that of the LSR. 

'For reference, the Milky Way's thick disk is described by (an, a^.az) = (63 ± 6, 39 it 4, 39 ± 4) km/s, iij^g = 51.5 km/s llSoubiran et al.l 
[2003) . and = 0.90 ± 0.18 kpc llJuric et al.ll2008f) . 



